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Abstract 

If large extra dimensions exist, the Planck scale may be as low as a TeV and mi¬ 
croscopic black holes may be produced in high-energy particle collisions at this energy 
scale. We simulate microscopic black hole formation at the Large Hadron Collider and 
compare the simulation results with recent experimental data by the Compact Muon 
Solenoid collaboration. The absence of observed black hole events in the experimental 
data allows us to set lower bounds on the Planck scale and various parameters related to 
microscopic black hole formation for a number (3 — 6) of extra dimensions. Our analysis 
sets lower bounds on the fundamental Planck scale ranging from 0.6 TeV to 4.8 TeV 
for black holes fully decaying into Standard Model particles and 0.3 TeV to 2.8 TeV for 
black holes settling down to a remnant, depending on the minimum allowed black hole 
mass at formation. Formation of black holes with mass less than 5.2 TeV to 6.5 TeV 
(SM decay) and 2.2 TeV to 3.4 TeV (remnant) is excluded at 95% C.L. Our analysis 
shows consistency with and difference from the CMS results. 


1 Introduction 

The Standard Model of Elementary Particles (SM) is one of the most successful hy¬ 
potheses in physics [1]. However, the SM fails to explain the hierarchy problem, i.e., the 
huge gap between the electroweak scale (Mew TeV) and the observed Planck scale 
(Mpi 10^® TeV). The ADD model [2] provides a way to solve the hierarchy problem 
by introducing a number n of large, compactified spatial dimensions (LEDs). Gravitons 
can propagate in the D{= n-|-4)-dimensional space-time bulk. SM particles are confined 
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to the 4-dimensional hrane. Assuming compactification on a torus with equal radii R, 
the observed Planck mass Mpi is related to the 4-dimensional fundamental Planck mass 
M* by Mpj = . If R is sufficiently large, the fundamental Planck mass M* 

may be as low as a few TeVs. 

If the ADD model is realized in nature, strong gravitational effects should manifest 
themselves in physical processes at the TeV scale. Gravitational phenomena at the TeV 
scale could include, for example, graviton and Kaluza-Klein mode production in particle 
scattering[3] and even microscopic black hole formation [4]. The Large Hadron Collider 
(LHC), operating at a center of mass energy of several TeVs, can be used to probe the 
appearance of these new physical processes and shed light on the existence of large extra 
dimensions [1, 5, 6, 7, 8]. 

To date, experimental results have not conhrmed the existence of large extra dimen¬ 
sions [9, 10, 11, 12, 13, 14, 15, 16]. These null results set upper bounds on the size of the 
large extra dimensions, or equivalently, lower bounds on the the fundamental Planck 
scale. Models with one or two large extra dimensions have been ruled out [17, 18], such 
as the Fermi Large Area Telescope [19]. Constraints on space-times with three large 
extra dimensions from astrophysical and cosmological experiments are generally very 
stringent, although they typically suffer from large systematic errors. The observation 
of Supernova SN1987A sets a lower limit on Md of 2.4 TeV for n = 3 [20], where the 
reduced Planck mass Md is related to M* by [5], 
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Neutron star-derived limits constrain Md to be larger than 76 TeV for n = 3 [21]. 
Collider experiments provide less stringent, albeit more accurate limits on Md for n > 3 
from non-observation of perturbative processes. These limits (in units of TeV) are shown 
in Table 1. 


Collider experiments 
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1.20 

1.15 

0.86 

0.80 

3.21 

1.25 

3.16 

0.98 

4.29 

1.20 

4.77 

4 

0.94 

1.04 

0.84 

0.73 

2.80 

1.26 

2.84 

0.96 

3.71 

1.17 

3.97 

5 

0.77 

0.98 

0.82 

0.66 

2.55 

1.26 

2.65 

0.92 

3.31 

1.12 

3.73 

6 

0.66 

0.94 

0.80 

0.65 

2.36 

1.29 

2.58 

0.88 

3.12 

1.07 

3.53 


Table 1: 95% C.L. lower limits on Md (TeV) from collider experiments. 


Lower bounds on the Planck scale can also be derived by non-observation of pro¬ 
duction and decay of TeV BHs in collider experiments and cosmic ray observations 
[34, 35, 36, 37, 38, 39, 40, 41]. The Compact Muon Solenoid (CMS) and A Toroidal 
LHC Apparatus (ATLAS) collaborations have conducted searches for BH signatures at 
the LHC [9, 42], setting limits on the production cross section and the minimal BH 
mass Mmin, i.e., the minimum mass at which a BH can form. Depending on model 
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assumptions, the CMS (ATLAS) collaboration excludes below 4.3 to 6.2 TeV (4.8 
to 6.2 TeV) at 95% C.L. 

The purpose of this paper is to revisit and extend the above results from CMS. 
We simulate production and decay of microscopic BHs at the LHC with the Monte 
Carlo generator CATFISH (v2.10) [18]. We then derive experimental bounds on the 
fundamental Planck mass and the number of extra dimensions by comparing simulation 
results with model-independent experimental limits on BH production from the CMS 
Collaboration [9]. The absence of observed black hole events in the CMS experimental 
data allows us to set bounds on various physical parameters of the ADD model and con¬ 
strain the minimum mass of TeV-scale black holes that may form in hadronic scattering 
processes. 


2 Black Hole Formation in Particle Collisions 


According to the Hoop Conjecture [43], a BH forms when a mass M is confined into a 
region of typical size equal to the Schwarzschild radius for that mass, i?s(M). Therefore, 
if two particles collide with center of mass energy ^/s and impact parameter smaller than 
Rs{y/s), a BH may form. If Rg <C R, as expected in the ADD scenario, the newly formed 
BH lives in a D-dimensional space-time with negligible curvature at the BH scale. In 
this case, the Schwarzschild radius of the BH can be expressed as [44, 45, 46], 


Rs = 
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where M = (1 — y)\/s and y is the fraction of energy which escapes into the bulk as 
gravitons, depending on the impact parameter. The Hoop Conjecture implies a BH 
production cross section cr{s,n,y) = ttFR'^, where the form factor D < 1 is related 
to y and accounts for the energy of the colliding particles that is not trapped in the 
event horizon, the so called “graviton energy loss at formation.” Since BH production 
in hadron colliders occurs at the parton level, the total cross section for a hadronic 
collision is obtained by integrating on the Parton Distribution Functions (PDFs) of the 
hadrons [45] 


>bh('S, u, y) — ^ ^ 
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X' 
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where fi(x,Q) are the PDFs with four-momentum transfer squared Q and z is the 
impact parameter normalized to its maximum value. The cutoff at small x is Xm = 
/{^[(l — y{z)\^}j where Mmin is the minimum-allowed mass of the BH. The total 
cross section in absence of graviton energy loss at formation is recovered by setting 
F = 1 (Black Disk (BD) cross section) [18]. 

If the initial BH mass is much larger than the Planck mass, a semiclassical treat¬ 
ment suggests that the newly-formed BH decays through four, possibly distinct stages: 
balding, spin-down, thermal evaporation and quantum decay [47]. During the balding 
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stage, the BH radiates multipole momenta and quantum numbers [6, 48], eventually 
settling down to a D-dimensional Kerr geometry. Angular momentum is radiated during 
the spin-down stage [48]. The Schwarzschild BH then decays into elementary particles 
through the Hawking mechanism (thermal evaporation stage). Most of the energy of the 
BH is radiated in this stage, with SM particles dominating the decay products. When 
the mass of the evaporating BH approaches the Planck scale, Qmin ~ Af*, the BH enters 
the quantum phase, where the decay ceases to be semiclassical and becomes dominated 
by quantum gravitational effects. 


3 BH Event Simulations 

Several Monte Carlo generators for BH production at colliders have been developed over 
the years: TRUENOIR [49], CHARYBDIS2 [50, 51], BlackMax [52, 53] (used by CMS), 
QBH [54] and an unnamed generator by Tanaka et al. [55]. Our analysis is based on 
CATFISH {Collider grAviTational Field Simulator for black Holes) [18]. CATFISH is 
a Fortran 77 Monte Carlo generator designed specifically for simulating BH events at 
CERN’s LHC. It incorporates three models for BH formation and cross section: BD, 
Yoshino-Nambu (YN) Trapped Surface (TS) [56], and Yoshino-Rychkov (YR) improved 
TS model [57]. The lack of a quantum theory of gravity requires a phenomenological 
treatment of the final stage. CATFISH offers the choice of simulating the quantum 
phase by either non-thermally decaying the BH into a number Up of hard quanta, each 
with energy Qmin/np, or forming a BH remnant. CATFISH also incorporates several 
other physical effects, such as inelasticity, exact field emissivities and corrections to 
semiclassical BH evaporation. The generator interfaces to the PYTHIA Monte Carlo 
fragmentation code [58] using the Les Houches interface [59]. In our analysis, we run 
CATFISH (v2.10) with the CTEQ6PDF PDF set and PYTHIA (v6.425) Tune Zl. 

The simulation of a BH event in CATFISH follows these steps. First, CATFISH 
computes the total and differential cross sections for the BH formation. The initial BH 
mass is sampled from the differential cross section. The BH is then decayed through 
the Hawking mechanism until the BH mass reaches the quantum limit, where a final 
non-thermal hard event is generated or a BH remnant is created. The unstable quanta 
emitted by the BH are instantaneously hadronized or decayed by PYTHIA, which also 
simulates initial- and final-state radiation particles. To determine the physics of BH 
formation and decay, CATFISH uses several external parameters and switches: 

• ADD parameters 

1. Fundamental Planck mass: ]yiSTAR=M*. 

2. Number of extra dimensions: NEXTRADIM= n, NEXTRADIM = 3, 4, 5, 6. 

• BH formation parameters 

1. Graviton energy loss at formation: GRAVITONLOSS = 0 (BD model), 1 (YR or 
YN TS models, see below). 

2. Gravitational loss model: GRAVITDNMODEL = 0 (YN TS model [56]), 1 (YR 
improved TS model [57]). 
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3. Minimum initial BH mass in Planck units: XMIN= Mmin/M*, XMIN> 1. 

• BH evaporation parameters 

1. Minimum quantum BH mass in Planck units: QMIN= Qmin/M^i QMIN < XMIN. 

2. Number of quanta emitted in the Planck phase: NP = Up. When NP=0, the 
BH forms a stable remnant with mass Qmin- 

The primary goal of this investigation is to determine lower bounds on the D- 
dimensional fundamental Planck scale for different values of NEXTRADIM, XMIN and NP 
using the model-independent 95% C.L. upper limits on the BH cross section from the 
CMS search [9]. We also derive lower limits on M m^n and XMIN for fixed MSTAR, NEXTRDIM 
and NP. For simplicity, we consider only BH formation with GRAVIT0NLDSS=0 (BD cross 
section), XMIN=QMIN and final decay in 2, 4, or 6 quanta or formation of a stable BH 
remnant (NP=0). The results with GRAVITONLOSS = 1 will be presented in a future 
report. 

The CMS search for BH events looks at excess transverse energy w.r.t. SM back¬ 
ground predictions [9]. The transverse energy St of an event is defined as the scalar 
sum of the transverse energies of all the final-state objects in excess of 50 GeV, i.e., jets, 
muons, electrons and photons satisfying the selection criteria discussed in Ref. [9]. The 
missing transverse energy is defined as the magnitude of the vector sum of the trans¬ 
verse momenta of all the final-state objects. If it is greater than 50 GeV, the missing 
transverse energy is added to St- The event multiplicity N is defined as the number of 
final-state objects which are used to calculate 5 t. BH events are expected to have high 
multiplicities. 


4 Results 

Lower bounds on M* and Mmin are derived by evaluating the partial cross section 
cr(5'T > 5™“) for events whose S't > and whose multiplicities are greater than 

some chosen value, given by, 

Cj(5t > = k ■ crpp_s.BH, (4) 

where S™™ is the minimal transverse energy chosen, and k is 

^ Num. of events with St > -S’™ 

Total Num. of events 

The behavior of the total cross section cr^p^BH as a function of M* and Mmin (he., XMIN) 
follows from Eqs.(2, 3). There are 3 factors: 1) The BD cross section cjbd is inversely 
proportional to a power of M*; 2) Since XMIN appears as a lower limit of integration 
in the total cross section Upp^-Qn, the greater XMIN, the smaller cjpp^BH is at a fixed 
M*; 3) In addition, the PDFs fall off rapidly at high Q. Taking into account all these 
factors, iTpp^BH is expected to decrease as M* (XMIN) increases at fixed XMIN (M*). The 
ratio k in Eq.(5) can be estimated by integrating over the spectra of visible final state 
particles over the range St > S™™ and then averaging over all final state particles, k is 
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an increasing function of the Hawking temperature, which, in turn, is a monotonically 
increasing function of M*. Thus, as M* increases, the graph of k vs. S'™™ flattens. 
Following the CMS collaboration [60], we choose the signal acceptance to be 100%. In 
summary, cr(ST > S™™) x H is expected to decrease as either M* or XMIN increases. 

The partial cross section (t(St > S™™) as a function of M* is obtained by running 
CATFISH with fixed NEXTRADIM, GRAVITONLOSS, XMIN = QMIN, NP. The results are shown 
by the upper two graphs in Figure (1) for NEXTRADIM = 3, XMIN = QMIN = 5, NP = 0, 
4. The lower two graphs in Figure (1) display it(St > S™™) as a function of XMIN for 
M* = 1.5 TeV, NEXTRADIM = 5, and NP = 0, 4. As expected, Figure (1) shows that the 
cross section decreases as M* and XMIN increase. 

Figure (1) can be used to determine bounds on M* and XMIN by comparing the 
simulated (t(S't > S™™) with the experimental limits (blue curves in the graphs). For 
example, the upper right plot shows that the lower limit on M* lies in the range 1.0 
TeV—1.2 TeV, and the bottom right plot shows that the lower limit on XMIN lies in the 
range 3.7—3.9. Multiplying this range by M* = 1.5 TeV, we obtain the lower limits on 
Mm\n = 5.55 TeV—5.85 TeV. 

We run CATFISH over a large range of parameter space as discussed in Appendix 
A. Figure ( 2 ) shows the exclusion region for M*. As expected, the lower limit M*^exp 
is a decreasing function of XMIN. The value of M*^exp does not strongly depend on NP, 
as long as NP % 0. If NP = 0, the bound on M*^exp becomes much smaller. This is due 
to the high transverse momentum of the BH remnant, which contributes to the missing 
energy. The lower limits on M* set upper bounds on Mmin (XMIN). As experimental data 
exclude values of M* < 1 TeV[18], our results for NP = 0 set an upper limit XMIN < 2.5. 
NP 7 ^ 0 results give the milder constraint, XMIN < 6 . More experimental limits are shown 
in Table 1, leading to more stringent constraints. For instance, CMS searches for events 
with an energetic jet and an imbalance in transverse momentum at y/s = 8 TeV [30] set 
the lower limit M* ~ 2.71 TeV for n = 3. The upper left panel in Figure (2) shows that 
the events with microscopic BHs decaying to remnants (NP = 0) are excluded, and the 
experimentally allowed range of XMIN is restricted to 1 ~ 2 for NP % 0. The lower limits 
of M* for n = 4,5, 6 do not exclude events with remnants as BH final products, but 
never the less, set strong constraints on the ranges where the semi-classical treatment 
is valid. 

Figure (3) shows the exclusion region for XMIN. As expected, the lower limit of XMIN 
is a decreasing function of M*. Figure (3) also shows that these limits do not depend 
on NP strongly when NP 7 ^ 0, but become much smaller at NP = 0. This figure can be 
combined with Figure (2) to constrain M* further. For example, if there are 3 extra 
dimensions, and a BH decays into 2 quanta in the quantum phase, the dashed curve 
(NP = 2, red onlne) in the upper left plot shows that XMIN > 4 at M* ~ 1.5 TeV. At 
the same time, the dashed curve in the upper left plot (n = 3) of Figure (2) indicates 
that if XMIN = 4, M* >1.5 TeV. This shows that the lower limits of M* are consistent 
with those of XMIN. In Figure (3), the ranges of M* were chosen in order to compare the 
CATFISH results with those of BlackMax and CHARYBDIS2 from Ref.’s [9, 42]. 

Figure (4) compares the lower limits of Mmin predicted by CATFISH with those 
from BlackMax and CHARYBDIS2 done by CMS Collaboration [9]. It shows that 
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as long as NP = 0, CATFISH’s limits are much smaller than those of BlackMax and 
CHARYBDIS2. When NP ^ 0, the situation is more complicated. CATFISH agrees 
with BlackMax (nonrotating BH model) very well when n = 4, but gives higher limits 
than BlackMax if n = 6. CATFISH gives limits similiar to those of CHARYBDIS2 
(nonrotating BH model) for n = 4, 6. While CHARYBDIS2’s stable remnant model 
with YR loss predicts smaller limits than CATFISH at n = 4,6, CHARYBDIS2’s boiling 
remnant model with YR loss gives smaller limits at n = 4, and similar limits to CATFISH 
at n = 6. Figure (5) shows the comparison between CATFISH predictions with those 
from ATLAS Collaboration [42]. CATFISH agrees with BlackMax and CHARYBDIS2 
when NP ^ 0, but predicts much smaller limits when NP = 0. The similarity displayed 
in the two figures is due to the fact that the 3 generators incorporate the same basic 
physics of microscopic BH formation and decay, but they differ from each other in the 
implementation of the quantum phase and/or taking graviton energy loss at formation 
into account. 


5 Conclusions 

In this work, lower limits on the fundamental Planck scale M* and minimal BH mass 
Mmin at formation have been obtained in a vast parameter space, using experimental 
upper limits on the partial production cross section of microscopic BHs [9]. Various 
models implemented in CATFISH have been explored and set different limits. In general, 
BH remnant models give milder constraints than non-remnant models. The predicted 
lower limits constrain the ADD model but do not exclude it. Future investigations will 
focus on performing a similar analysis to that carried out in the present paper with the 
additional feature of graviton energy loss during BH formation (GRAVITONLOSS = 1). 
Another refinement of the models is to include the effects of the generalized uncertainty 
principle [61]. The steps outlined above for the simulation of microscopic BH events can 
also be carried out for string balls, string resonances and other exotic particles. 
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A Appendix A: Procedure for determining 
and Mmin limits 


The primary goal of the work is to determine the limits on the L)-dimensional funda¬ 
mental Planck scale M*, or equivalently, the reduced Planck scale Md used in Ref. [9], 
for different NEXTRADIM (=3, 4, 5, 6), GRAVITONLOSS (=0), XMIN (= QMIN) and NP (=0, 
2, 4, 6). The model-independent 95% C.L. upper cross section limits for counting ex¬ 
periments (Figures 6 and 7 in Ref. [9]) were used. CATFISH accepts only XMIN > 1, so 
the lower limit of XMIN is 1 and its upper limit is determined by noticing that the mass 
of the BH should be less than the center of mass energy of the LHC, 

XMIN X MSTAR < 8 TeV. (6) 


The stepsize of XMIN is set to 0.5. We use a bisection method to find the limits on the 
D-dimensional fundamental Planck scale: 

Step 1. The possible range of M* is determined by fixing all other parameters. The 
lower limit is mo = 1 TeV and the upper limit, Mq, satisfies Eq. (6) and the condition 


a ■ L > 


N, 


run 


(7) 


where a is the cross section of the production, L = 12 fb ^ is the LHC integrated 
luminosity, and Wun is the number of events of each run. In our simulations we choose 
Vrun = 10^. 

Step 2. CATFISH is run with MSTAR=Mi = (mo -|-Mo)/2 and the simulation output 
is compared with CMS results [9]. 

Step 3. If Ml is allowed by experimental data, i.e., the simulated cross section (times 
the detector acceptance A) for BH production is too small, M* must be smaller than 
Ml. CATFISH is run with MSTAR=M 2 = (mo -|- Mi)/2. Otherwise, CATFISH is run 
with MSTAR=M 2 = (Ml -k Mo)/2. 

Step 4. Step 3 is repeated i times until |Mj — Mj_i| < AM, where Mi is the result 
of the i-th simulation and AM is the required precision, AM = 0.1 TeV. 

Step 5. The more likely value of M* is then determined as 


Mjf ^ exp 


Mi + Mj_i 

2 


The error on M*^exp is AM*^exp = p|Mj — Mj_i|, where p is determined at a given 
confidential level (C.L.) for the interval [Min(Mj, Mj_i),Max(Mj, Mj_i)], 

1 - C.L./lOO 
- 2 -• 

The lower limit on the Planck mass, M*^ expj is determined for a given choice of NEXTRADIM, 
GRAVITONLOSS, XMIN (= QMIN) and NP. Different sets of these parameters are chosen 
and Steps 1-5 are repeated to determine M*_exp as function of the parameters. Simi- 
larily, the lower limit of Mmin or XMIN can be obtained for different choices of MSTAR, 
NEXTRADIM, GRAVITONLOSS (=0) and NP. 
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Figure 1: > S'™”) x A as a function of M* (Upper two graphs, numbers being M*’s chosen in 

units of TeV) or XMIN (Lower two graphs, numbers being XMIN’s) at NP = 0 (Right two graphs) 
and NP = 4 (Left two graphs). The model-independent 95% CL experimental upper limits for 
counting experiments from CMS Collaboration are also shown. The multiplicity is > 3. 
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Figure 2: Simulated lower limit on M* vs. XMIN as a function of NP and NEXTRADIM = 3 
left), 4(top right), 5 (bottom left), and 6 (bottom right). 
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Figure 3: Simulated lower limit on XMIN vs. M* as a function of NP and NEXTRADIM = 3 
left), 4(top right), 5 (bottom left), and 6 (bottom right). 
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Figure 4: Comparison of the predictions on lower limits of from CATFISH with those from 

BlackMax and CHARYBDIS2 at NEXTRADIM = 4 (left) and 6 (right). The results of BlackMax 
and CHARBDIS are extracted from Fig. 4 in Ref. [9]. 
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Figure 5: Comparison of the predictions on lower limits of from CATFISH with those from 

BlackMax and CHARYBDIS2 at NEXTRADIM = 4 (left) and 6 (right). The results of BlackMax 
and CHARBDIS are extracted from Fig.’s 8 and 10 in Ref. [42]. 
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